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The detrending moving average (DMA) algorithm is one of the best performing methods to 
quantify the long-term correlations in nonstationary time series. Many long-term correlated time 
series in real systems contain various trends. We investigate the effects of polynomial trends on 
the scaling behaviors and the performances of three widely used DMA methods including backward 
algorithm (BDMA), centered algorithm (CDMA) and forward algorithm (FDMA). We derive a 
general framework for polynomial trends and obtain analytical results for constant shifts and linear 
trends. We hnd that the behavior of the CDMA method is not influenced by constant shifts. In 
contrast, linear trends cause a crossover in the CDMA fluctuation functions. We also hnd that 
constant shifts and linear trends cause crossovers in the huctuation functions obtained from the 
BDMA and FDMA methods. When a crossover exists, the scaling behavior at small scales comes 
from the intrinsic time series while that at large scales is dominated by the constant shifts or linear 
trends. We also derive analytically the expressions of crossover scales and show that the crossover 
scale depends on the strength of the polynomial trend, the Hurst index, and in some cases (linear 
trends for BDMA and FDMA) the length of the time series. In all cases, the BDMA and the FDMA 
behave almost the same under the influence of constant shifts or linear trends. Extensive numerical 
experiments confirm excellently the analytical derivations. We conclude that the CDMA method 
outperforms the BDMA and FDMA methods in the presence of polynomial trends. 

Keywords-. Fractal Analysis; Detrending Moving Average (DMA); Scaling law; Crossover Behav¬ 
ior; Polynomial Trend; Constant Shift; Linear Trend. 


wxzhou@ecust.edu.cn 



2 


I. INTRODUCTION 


Many natural, social and technological systems exhibit complex behavior characterized by long-term power-law 
correlations There are a wealth of methods developed to determine the correlation strength in lo ng -term correlated 
time series [JQ- The most classic method is Hurst analysis or rescaled range analysis (R/S) 0, 0 ■ Other popular 
methods include wavelet transform module maxima fWTMM) approaches 0, lol-H^ , detrended fluctuation analysis 
(DFA) 0 based on the fluctuation analysis (FA) Ml: detrending moving average analysis (DMA) [TMl based 
on the moving average or mobile average techniqu e |19l| . and so on. These methods have been generalized in many 
directions, such as objects in high dimensions 01117 detrended cross-correlation analysis and its variants for two 
time analysis JiQ: detrended partial cross-correlation analysis for multivariate time series [3ll - [33l| . multifractal 
analysis [3J, 0 , and so on. 

An important issue is to compare the performance and relative merits of these estimators, which has been conducted 
through extensive numerical experiments. With time series generated from the modified Fourier filtering method [s^ . 
Xu et al. found that DFA is superior to different DMA variants [s^- Bashan et al. observed that CDMA performs 
comparably well as DFA for long time series with weak trends and slightly outperforms DFA for short data with weak 
trends [s^. Based on fractional Gaussian noises (FGNs) generated from the Davies-Harte algorithm and fractional 
Brownian motions by summing the FGNs, Serinaldi found that DFA and DMA have comparable performances [40| . 
Jiang and Zhou reported that DFA and CDMA perform similarly and both of them outperform the BDMA and 
FDMA methods [1^, in which the FBMs are generated using the Fourier-based Wood-Chan algorithm [4l|. Huang 
et al. reported comparative performances of FA and DFA for FBMs with = 1/3 [i^, which were generated with 
the Wood-Chan algorithm [^. Bryce and Sprague reported that FA outperforms DFA, for FGNs with H = 0.3 
[ 4 ^ . which were generated using the Davies-Harte algorithm [s^, while Shao et al. found that CDMA has the best 
performance, DFA is only slightly worse in some situations, and FA performs the worst [3. It is not unreasonable 
that the conclusions are mixed because different studies used different time series generators and different lengths. 

Time series in real complex systems usually contains various forms of trends and nonstationarity. Hence, another 
important issue concerns the effects of trends and nonstationarity on the scaling behaviors of different methods. 
Montanari et al. investigated the effects of periodicity on several methods such as aggregated variance method, 
Higuchi’s method, R/S analysis, periodogram method. Whittle method, and so on 0]. Kantelhardt et al. studied 
the effects of polynomial trends and oscillatory trends on the different orders of DFA [d^. Hu et al systematically 
studied the effects of linear, periodic, and power-law trends on DFA [i^. Chen et al. considered the presence of non- 
stationarity and nonlinear filters in the DFA analysis [13, researched the effect of missing extreme data 

on DFA [ 4 ^ . Song and Shang investigated the effects of five trends based on linear and nonlinear filters on multifractal 
DCCA based on DFA [ 5 ^. In most cases, a crossover appears in the scaling behavior of the DFA fluctuation functions, 
which makes it difficult to estimate the intrinsic long-term correlations in time series. Many efforts have been made 
to reduce or minimize these effects on the DFA method 

However, studies on the effect of trends on the detrending moving average analysis are rare, although DMA is “The 
Method of Choice” as DFA [IJ. To our knowledge, one such study is to minimize the effect of period trends on 
the DMA method [s^]- In this work, we aim at contributing this direction by investigating the effects of polynomial 
trends on the scaling behavior of DMA methods. We derive analytically the results for constant shift and linear trend 
and confirm these results using numerical experiments. 


II. DETRENDING MOVING AVERAGE ALGORITHMS 


The algorithms of the detrending moving average analysis are described as follows (l6l - [l^ 1^ . [591 - 0 . 

Step 1. Consider a time series x{t), t = 1, 2, • • • ,N. We construct the sequence of cumulative sums 

t 

X{t) = Y,x{i), t = (1) 

i=l 


Step 2. Consider a box [t — si,t + S2] of size s = si -I- S2 -f 1, where si = [(s — 1)(1 — 0)], S2 = [(s — 1)0J, L^J is 
the largest integer smaller than x, [a;] is the smallest integer larger than x, and 6 is the position parameter with the 
value varying in the range [0,1]. Calculate the moving average function X{t) in a moving window El: 

X{t) = - f2 Xit-k). 

k= — S2 


( 2 ) 
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Hence, the moving average function considers si data points in the past and S 2 points in the future. We consider 
three special cases in this paper. The first case 9 = 0 refers to the backward moving average [13 , in which the 
moving average function X{t) is calculated over all the past n — 1 data points of the signal. The second case 9 = 0.5 
corresponds to the centered moving average [s^ - where X{t) contains half past and half future information in each 
window. The third case 0 = 1 is called the forward moving average, where X{t) considers the trend of n — 1 data 
points in the future. 

Step 3. Detrend the signal series by removing the moving average function X(i) from X(i), and obtain the residual 
sequence e(i) through 


e(z) = X(t)-X(2), (3) 

where n — [(s — 1)0J ^ i ^ X — [(s — 1)0J. 

Step 4- The residual series e(i) is divided into W disjoint segments with the same size s, where Xg = [^X/n — IJ. 
Each segment can be denoted by e„ such that ev{i) = + *) for 1 ^ i ^ s, where I = (v — l)s. The root-mean-square 

function Fy{s) with the window size s can be calculated by 

FS{s) = -j24i^)- ( 4 ) 

Step 5. Varying the values of segment size s, we can determine the power-law relation between the function F(s) 
and the size scale s, which reads 


F{s) - bs^. 


( 5 ) 


III. POLYNOMIAL TRENDS 


Consider a signal composed of a signal x(t) with zero mean and an additive trend u{t) 

z{t) = x{t) + u{t) 

The profile of z{t) is the sum of the profiles of x{t) and u{t): 

Z{t) = Xit) + U{t) 

and the moving average at time t is 

Zit) = X{t) + U{t) 

When q = 2, the overall fluctuation is 

N 

F^{s)=Y,[z{t)-z{t)r 

N 

= Y^[xit)-x{t)+u{t)-u{t)r 

N 

= T'x(s) + F^{s) + 2y^jex{t)eu{t)] 

i=l 

where ex{t) = X{t) — X{t) and eu{t) = U{t) — U{t). If ex{t) and e„(t) are uncorrelated, we have 

F^As) = Fl{s) + Fl[s) 

which is the superposition rule (d^ . 

We consider polynomial trends added to the increments series: 

m 

u{t) = ^ apfP 

p—0 


( 6 ) 

( 7 ) 

( 8 ) 


(9) 


( 10 ) 


( 11 ) 
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and the cumulative sum is 


U{t)=Y^ u{i) = ^ap^tP 

i—1 p—0 


( 12 ) 


According to Faulhaber’s formula, the sum of powers X)i=i expressed as follows: 




(13) 


k=0 


where the coefficients B^. are the Bernoulli numbers. For p=0, 1, 2 and 3, we have X)i=i X)i=i + ^)/2i 

J2l=i + t )/6 and X)i=i + f^)/4- It follows that 


U{t) — ( tto + -ai + -02 ) ^ + ( + 7 ; 0,2 + - 70,3 ) t 
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+ ( ^02 + -03 ) + - ast ^ 


(14) 


The moving average at t from t — si to t + S 2 is 


t + S2 


z{t) = - J2 z{t-k) = - Y. z{k) 


(15) 


k— — S2 


k—t — si 


where s = si + S 2 — 1 is the window size. 


IV. CONSTANT SHIFT: THE CASE OF p = 0 
A. Analytical results 


In this case, we consider oi = 02 = 03 = 0. The trend is a constant shift 

u{t) = ao 

The cumulative sum, or the profile, is 


ui.t) = y^opt 


2=1 

The moving average at t obtained from t — si to t S 2 for window size s is 

U{t) =ao(^t+ - Si 

where si + S 2 + 1 = s. Since Si = (s — 1)(1 — 6 ) when 0 = 0,0.5 and 1 (note that s should be odd), we have 

(20-l)(s-l)- 


U{t) = ao 


t + ■ 


and the residual of the trend after removing the moving average is 

(20-l)(s-l) 
e„(t) = oo- - - 

which is a constant for a given window size s. Hence the superposition rule holds. 
When 9 = 0.5, we have 


(16) 

(17) 

(18) 

(19) 

( 20 ) 


FHs) 


1 

N 


N 




( 21 ) 
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FIG. 1. Effect of constant shift on the DMA algorithms. Each curve in (a-d) represents a fluctuation function averaged over 
50 repeated simulations, (a) Log-log plots of {F) of the FGNs with different Hurst indexes and the FGNs with constant shift 
ao — 0.2 in the increments against s for the CDMA method. The Hurst index H varies from 0.1 (bottom) to 0.9 (top) with 
a step of 0.1. (b) Log-log plots of (F) against s using the BDMA and FDMA methods for iLin = 0.1 and ao = 0.0012. (c) 
Log-log plots of {F) against s using the BDMA and FDMA methods for H — 0.5 and ao — 0.0176. (d) Log-log plots of {F) 
against s using the BDMA and FDMA methods for H = 0.9 and ao — 0.257. (e) Power-law dependence of the crossover scale 
Sx on the constant shift ao using the BDMA method for different Hurst indexes, varying from 0.1 (left) to 0.9 (right) with 
a step of 0.1. (f) Validation of Eq. (1241) . The crossover exponent a is the power-law exponent in (e) and K = —1/(1 — H) 
and k = —1/(1 — h), where FI is the input Hurst index for the generation of FGNs and h is the estimated Hurst index of the 
generated FGNs using BDMA. 


which is independent of the constant shift term oq. It indicates that, if a; is a fractional Gaussian noise, there is no 
crossover in the scaling. 

When 0 = 0 and 0=1, the detrended fluctuation is 

+ ( 22 ) 

which depends on ao and s. There is a crossover s = Sx in Fz{s). For s < Sx, the behavior of F/(s) is very close to 
the behavior of i/c(s), while for s > Sx, the behavior of ^/(s) is very close to the behavior of F„(s). The crossover 
scale Sx is the solution to the following equation 


Fx{s) = Fu{s). 


It follows that bs^ = ao(s — l)/2. When s 1, we have 


Sx 




which shows that Sx is a power-law function of oq with the exponent being —1/(1 — H). 


(23) 


(24) 


B. Numerical experiments 

We perform numerical simulations to verify the correctness of the main results, Eq. (|^ and Eq. (l24|) . derived in 
Sec. IIV Al We employ the Davies-Harte algorithm [s^ to generate fractional Gaussian noise (EGN) with given Hurst 
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indexes H. There are other comparable generators such as the wavelet-based FMB generator and the random 
midpoint displacement algorithm [^. However, the Davies-Harte algorithm performs slightly better [d^. In our 
simulations, we consider different Hurst indexes H, which range from 0.1 to 0.9 with a step of 0.1. For each H, we 
generate 50 FGN time series x(t) of length 10®. The constant shift oq is added to each point of the each FGN series. 
The DMA fluctuation functions presented below are averaged over the 50 realizations. 

Figure [IK illustrates the averaged fluctuation functions obtained from the CMDA method for different H values. 
The two curves for the original FGN time series and for the shifted FGNs with qq = 0.2 overlap excellently. In 
addition, all the curves have excellent power-law forms with the slopes being the corresponding H values. Ghanging 
the value of Oq has no impact on the results. Therefore, Fig. UK verifies Eq. (1^ exactly. 

Figure [IK shows the average fluctuation functions of the constantly shifted FGNs with oq = 0.0012 and H = 0.1 
obtained from the BDMA method and the FDMA method. The two curves overlap nicely. The fluctuation functions 
exhibit a clear crossover. When s <C Sx, the fluctuation functions overlap with the fluctuation function of FGNs 
with the slope being H = 0.1. When s ^ Sx, the fluctuation functions overlap with the fluctuation function of the 
constant oq with the slope being H = 1. These observations are consistent with Eq. (1^^ . Eigure [IK and Fig. [TK 
show the results for H = 0.5 and oq = 0.0176 and for H = 0.9 and oq = 0.257, respectively. These results are also 
consistent with the prediction of Eq. 

We can determine the crossover scale Sx by two methods. The hrst one is to use Fj;{sx) = Fu{sx), determining 
the intersection point {sx, Fj;{sx)) of the solid line and the dashed line in each plot (Fig. [IK-Fig. [Ui). However, this 
method uses a priori information about the underlying FGNs and the constant shift. An alternative method is as 
described below. We pinpoint the point Sm on the fluctuation curve that is the farthest from the line connect the two 
endpoints. We perform a linear ht the first few point from the right point to a point that is in the middle of the right 
point and {sm, Fz{sm)) to obtain a first straight line and similarly a second straight line based on the right part of 
the Fz curve. The crossover scale is determined by the intersection of these two straight lines. In this procedure, the 
choices of the right point Sm for the left part of the fluctuation function and the left point of the right part of the 
fluctuation function can vary, which does not influence the determination of Sx- In our analysis, we simply use the 
five left-most data points and the hve right-most data points in the linear regressions and obtain the intersection of 
the two regressed lines treating as Sx- Note that there are 60 points in each fluctuation function. Figure [IK shows 
the dependence of Sx as a function of oq for different F[ values. For every F[ value, we observe a nice power-law 
relationship: 


Sx a^, (25) 

which is consistent with the power-law form expressed in Eq. (12411 . 

Eigure[If shows that the lines become steeper for larger Hurst indexes. We fit the data points for each F[ to estimate 
the power-law exponent a. We then define and calculate the following two quantities: 

K = H) (26) 

and 

k = -l/{l-h), (27) 

where F[ is the input Hurst indexes for the synthesis of the FGNs and h is the output Hurst indexes of the synthesized 

FGNs using the BDMA method. We plot K against a and k against a in Fig. [If. We observe that 

K^k^a, (28) 

except for FI = 0.9. Equations (1^51) and (1^51) verify excellently Eq. (IMl) . 

The choice of ag values are not arbitrary. Due to the finite size of the generated EGNs, too large oq will result in 
very small Sx so that the resulting fluctuation function becomes a straight line with the slope being 1, while too small 
ag will result in very large Sx so that the resulting fluctuation function becomes a straight line with the slope being 
i7. In both cases, the crossover cannot be recognized. In the numerical experiments, for each F[, we use 10 ag values 
that are evenly spaced in the logarithmic scale. For instance, ag values are distributed in [0.00018, 0.012] for H = 0.1, 
in [0.0064, 0.0621] for F[ = 0.5, and in [0.203, 0.346] for F[ = 0.9. In this way, the crossovers can be identified. The ag 
values used in Fig. [TK to Fig. [Ui are the fifth in each of the 10 ag values. 

It is clear that the numerical results illustrated in Fig. [U verify the analytical results in the previous subsection. 
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V. LINEAR TREND: THE CASE OF p = 1 
A. Analytical results 

We now consider the case of linear trends with ap = 02 = 03 = 0. The linear trend is 

u{t) = ait. 

The profile of u{t) is 

U{t) = ai{t^ + t)/2, 

and the moving average is 

~, /I 9 s — 2 si ssi s? 1 

where si = (s — 1)(1 — 9). When 0 = 0,0.5,1 (note that s should be odd), we have 


U[t) = ai 


1 2 2(1-0) + s(20-1) ' 

2^ ^ 2 


+ L, 


where 


L = ai(s - l)(s - 3s9 + 3s0^ - 30^ + 60 - 2)/6. 
eu{t) = U(t)- U{t) = At- L, 

A = -ai(s-l)(20-l)/2. 

1 ^ 

F^s) = Yl[<^x{t) + e„(t)]2 

Applying the superposition rule and Faulhaber’s formula, we have 

1 N 


The residual series is 


where 


The detrended fluctuation is 


F^{s) = j;^Yl^el{t) + el{t)] 


t=l 


= F^ + lf2{At-L) 


t = l 

2 2 A2(27V2+ 3A+1) , , 

= + L^ + —^^ - AL{N + 1) 


(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


(37) 


When 0 = 0.5, we have A = 0 and L = ai(s^ — l)/24. Inserting them into Eq. (|37)l . it follows immediately that 


f! = Fl 


ai(s^ - 1) 

24 


We obtain the crossover scale Sx as follow 


ai J 




(38) 


(39) 


We notice that Sx depends on m but not on N. 


















FIG. 2. Effect of linear trend on the CDMA algorithm. Each curve in (a-d) represents a fluctuation function averaged over 50 
repeated simulations, (a) Log-log plots of {F) of the FGNs with H = 0.1 contaminated by a linear trend with ai = 9.55 x 10“® 
in the increments, (b) Log-log plots of {F) of the FGNs for H = 0.3 and ai = 2.70 x 10~®. (c) Log-log plots of {F) of the FGNs 
for F[ = 0.7 and ai — 2.18 x 10“'^. (d) Log-log plots of (F) of the FGNs for FI = 0.9 and oi = 4.61 x 10“"^. (e) Power-law 
dependence of the crossover scale Sx on the linear coefficient ai for different Hurst indexes, varying from 0.1 (left) to 0.9 (right) 
with a step of 0.1. (f) Validation of Eq. (I39II . The crossover exponent a is the power-law exponent in (e) and K = —1/(2 — H) 
and k = —1/{2 — h), where H is the input Hurst index for the generation of FGNs and h is the estimated Hurst index of the 
generated FGNs using CDMA. 


When 0 = 0, we have L = ai(s — l)(s — 2)/6 and A = ai{s — l)/2. The detrended fluctuation is 


When 1 <C s -C Ai, we have 


There is a crossover at scale 


+ 


1 2 


ai(s — l)(s — 2) 

6 

al{2N'^ + m + l){s - If 
24 

a2(iV + l)(s-l)2(s-2) 

12 




a'lN'^s'^ 


12 


Ns N^\ 

12^12) 


Sx 





(40) 


(41) 


(42) 


which depends on ai and N. 
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When 0 = 1, we have A = —ai(s — l)/2 and L = ai(s^ — l)/6. It follows that 


The crossover Sx is derived as follows 


+ aM 


36 12 


12 




12 


Sx — 


( VT2b\ 

\a,N J 




(43) 


(44) 


which depends on ai and N. We notice that the crossover scales for the FDMA and the BDMA have approximately 
the same expression. 


B. Numerical simulations 

We now perform numerical simulations to verify the correctness of the main results derived in Sec. IV Al in particular 
Eq. (1391) for the CDMA method (0 = 0.5), Eq. (l42l) for the BDMA method (0 = 0) and Eq. (l44t for the FDMA method 
(0 = 1). The procedures of numerical experiments are the same as for the case of constant shift in Sec. IIVBI 

The results for the CDMA method are illustrated in Fig.j^J In Fig.[2Ji, we show the averaged fluctuation function of 
the FGNs with JI = 0.1 contaminated by a linear trend with ai = 9.55 x 10“®. We observe an evident crossover Sx in 
the fluctuation function. When s s x, the fluctuation function overlaps excellently with the fluctuation function of 
FGNs with the slope being H = 0.1. When s Sx, the fluctuation function overlaps excellently with the fluctuation 
function of the linear trend ait with the slope being H = 2. These observations are consistent with Eq. (l38)l . We 
present respectively the results for H = 0.3 and oi = 2.70 x 10“^ in Fig. [2b, for FI = 0.7 and ai = 2.18 x 10“^ in 
Fig. [2}:, and for H = 0.9 and ai = 4.61 x 10“^ in Fig. |2}1. All these results are also consistent with the prediction of 
Eq. (1551) . 

We adopt the same procedure as the case of constant shift in the determination of the crossover scale s x for different 
H values. For each FI, we choose 10 values of ai, which are evenly spaced in logarithmic scales. Figure |2b shows 
the dependence of Sx as a function of ai for different F[ values. For fixed ai, Sx increases with H, suggesting that 
stronger long-term correlation in the FGNs corresponds to wider scaling range in the intrinsic fluctuation function. 
For fixed H, Sx decreases with ai, indicating that stronger trend will narrow the scaling range of the intrinsic FGNs 
and make it more difficult to determine the intrinsic Hurst index. For every FI value, we observe a nice power-law 
relationship: 


Sx ^ a“, (45) 

which is consistent with the power-law form expressed in Eq. (1391) . 

Figure|2f shows that the lines become steeper for larger Hurst indexes. We fit the data points for each FI to estimate 
the power-law exponent a. We then define and calculate the following two quantities: 

K = -l/{2- H) (46) 

and 

k = -l/{2-h), (47) 

where FI is the input Hurst indexes for the synthesis of the FGNs and h is the output Hurst indexes of the synthesized 
FGNs using the CDMA method. We plot K against a and k against a in Fig. [Tf. We observe that all the points fall 
on the diagonal line 


K K, k K, a. 


(48) 


Equations (1451) and (l48l) verify excellently Eq. (l39|) . 

The results for the BDMA and FDMA methods are depicted in Fig.jS) Applying the BDMA and FDMA methods, 
we show in Fig.|3b the averaged fluctuation functions of the FGNs with 77 = 0.1 contaminated by a linear trend with 
oi = 2.28 X 10“®. The two curves overlap nicely. We observe an evident crossover Sx in the fluctuation functions. 
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FIG. 3. Effect of linear trend on the BDMA and FDMA algorithms. Each curve in (a-e) represents the fluctuation function 
averaged over 50 repeated simulations, (a) Log-log plots of (F) of the FGNs with H — 0.1 contaminated by a linear trend 
with ai = 2.28 x 10“® in the increments, (b) Log-log plots of (F) against s for the case ol H = 0.3 and ai = 7.53 x 10“®. (c) 
Log-log plots of (F) against s for the case of 77 = 0.7 and ai = 3.95 x 10“®. (d) Log-log plots of (F) against s for the case of 
77 = 0.7 and ai = 1.21 x 10“^. (e) Log-log plots of (F) against for the case of 77 = 0.9 and ai = 4.45 x 10~^. (f) Power-law 
dependence of the crossover scale Sx obtained from the BDMA method on the linear coefficient ai for different Hurst indexes, 
varying from 0.1 (left) to 0.9 (right) with a step of 0.1. (g) Validation of Sx ~ in Eq. (1421) . The crossover exponent a 

is the power-law exponent in (f) and K — —1/(1 — 77) and k = —1/(1 — h), where 77 is the input Hurst index for the generation 
of FGNs and h is the estimated Hurst index of the generated FGNs using BDMA. (h) Power-law dependence of the crossover 
scale Sx obtained from the BDMA method on the time series length N for different Hurst indexes, varying from 0.1 (bottom) 
to 0.9 (top) with a step of 0.1. (i) Validation of Sx ~ Eq, (|42|) . The crossover exponent /3 is the power-law 

exponent in (h) and K = —1/(1 — 77) and k = —1/(1 — h). 


When s <C Sx, the fluctuation function overlaps excellently with the fluctuation function of FGNs with the Hurst 
index H = 0.1, which is shown as a solid straight line. When s Sx, the fluctuation function overlaps excellently 
with the fluctuation function of the linear trend ait with the slope being 77 = 1, which is shown as a dashed straight 
line. We present respectively the results for H = 0.3 and ai = 7.53 x 10“® in Fig. [31 d, H = 0.7 and oi = 3.95 x 10“® 
in Fig. [3}:, for H = 0.7 and ai = 1.21 x 10“^ in Fig.[3Ji, and for H = 0.9 and Oi = 4.45 x 10“^ in Fig. [5J;. All these 
results are consistent with the prediction of Eq. dni) and Eq. (PI) . which have the same expression. Note that these 
two expressions are approximations of Eq. (1371) . The difference between the functions of the BDMA and FDMA 
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methods is about 2AL. These results show that this difference is ignorable and the approximation leading to Eq. dSD 
and Eq. (H51) is reasonable. 

We determine the crossover scale Sx for different H and oi values. Eigure [3f shows the dependence of Sx as a 
function of oi for different H values. It is again found that Sx increases with H for fixed oi and decreases with oi 
for fixed H. For every H value, we observe a nice power-law relationship: 

Sx ~ a?. (49) 

We determine the power-law exponents a for different H values and plot K = —1/(1 —i/) against a and k = —1/(1 —/i) 
against a in Fig. [5^. We observe that all the points fall on the diagonal line 

(50) 

except for H = 0.9. Equations (HHl) and (1501) verify excellently Eq. (|00|) . 

We now investigate the dependence of the crossover scale Sx on the length N of time series. In Fig. [3jr, we plot Sx 
against N in log-log scales for different Hurst indexes H and ai values: ai = 7 x 10“® for H = 0.1, ai = 2 x 10“^ 
for H = 0.2, ai = 4 X 10"^ for H = 0.3, ai = 6 x 10"^ for H = 0.4, ai = 1 x lO'^ for H = 0.5, ai = 2 x 10"® for 
H = 0.6, oi = 3 X 10“® for H = 0.7, oi = 5.4 x 10“® for H = 0.8, ai = 6.4 x 10“® for H = 0.9. In our numerical 
experiments, N ranges from 50000 to 100000 with a step of 5000. The determination of this range is not arbitrary. 
If we include shorter time series, say N ~ 10^, we have Sx 10® ®, which is larger than 7V/10 and thus cannot be 
detected in the fluctuation function. If we include longer time series, say N ^ 10®, we have Sx 10^"^^, which is 
again hard to identify. 

For every H value, we observe a nice power-law relationship: 

(51) 

The power-law scaling for H = 0.9 is the worst. When N is small (left part) or large (right part), the crossover 
scale is large or small, which makes it very difficult to identify because the crossovers are near the end points of the 
fluctuation function. 

We determine the power-law exponents /3 for different H values and plot K = —1/(1 — 77) against j3 and k = 
— 1/(1 — h) against /3 in Fig.[5J. We observe that all the points fall on the diagonal line 

KPikPi (3, (52) 

except for H = 0.9. Equations (1511) and (15^ verify excellently Eq. (H^ . 


VI. SUMMARY 

In this paper, using fractional Gaussian noises (FGNs) with different Hurst indexes, we have investigated the effects 
of polynomial trends on the scaling behaviors and the performance of three widely used DMA methods including 
backward algorithm (BDMA), centered algorithm (GDMA) and forward algorithm (FDMA). We derived a general 
framework for polynomial trends and obtained analytical results for constant shifts and linear trends in the FGNs. 
We performed extensive numerical experiments which confirm excellently the analytical derivations. 

We first considered constant shifts in the FGNs. We found that the behavior of the GDMA method is not influenced 
by constant shifts. In contrast, constant shifts result in a crossover Sx in the fluctuation function when the BDMA 
and FDMA methods are applied. The crossover Sx scales as a power law of the constant shift, which increases with 
the Hurst index H and decreases with the strength of constant shift oq. 

We then considered linear trends in the FGNs. We found that a linear trend ait causes a crossover in the fluctuation 
function of the CDMA method. The crossover scale Sx is a power law of the strength oi of the linear trend, which 
increases with the Hurst index H and decreases with the strength of linear trend oi. A linear trend also results in a 
crossover in the fluctuation function when the BDMA and FDMA methods are applied. The crossover Sx scales as 
a power law of the production of the strength of linear trend and the length of the time series, which increases with 
the Hurst index H and decreases with the strength of constant shift ai and the length of linear trend. It is intriguing 
that longer time series are less resistent to linear trends when the BDMA and FDMA methods are adopted. 

When a crossover appears, the left part of the fluctuation function with the scales less than Sx reflects the behavior 
of FGNs, while the right part with the scales greater than Sx is dominated by the polynomial trend. Our findings 
show that time series with larger Hurst indexes are more resistent to polynomial trends and the polynomial trend 
has the same effect on the BDMA and the FDMA. Because a large crossover scale will make it easier to estimate 
the intrinsic Hurst index, we conclude that the GDMA method outperforms the BDMA and FDMA methods in the 
presence of polynomial trends (see Fig.[T^ versus Fig.[Tj3 for constant shifts and Fig.[2}3 versus Fig. [3]; for linear trends). 
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